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1. Introduction. 

The problem of the origin of the cosmic rays is still uncer- 
tain. Since the observed cosmic ray composition, spectrum 
and anisotropy involve a lot of factor, for example, the 
spectrum of the primordial cosmic rays, the mechanism of 
acceleration and propagation as well as the distribution of 
the interstellar medium and magnetic field. Unfortunately, 
all of these are remained open, though people gets more and 
more information and ideas. As a theory, it should explain 
the support of particles and energy, the mechanism of acce- 
leration and propagation as well as some important features 
obtained directly from cosmic ray experiments, such as the 
power spectrum and the knee at about lO^eV 1 ' 1 , the near cons- 
tant amplitude of the first harmonic of anisotropy among 
10"-10“*eV and the amplitude varying as E°* s above ExlO'^eV^,^ 
And it should also account for the relative abundance and 
the flux of cosmic rays. But so far there is no model which 
can interpret all of these phenomena. , 0 

In a general opinion the cosmic rays of 10* -10 eU are 
the galactic origin, but above they are the extragalactic 
origin. People also acceptes that the cosmic rays propagate 
in a diffusion way, variaty of diffusion mechanism derive 
different diffusion coefficient as a function of energy. 

But anyhow diffusion always makes the energy spectrum 
steeper. 

There are two kinds of models for interpreting the tlJ 
knee of the cosmic ray spectrum. One is the leaky box model 
in which every nuclei has different escape starting ri- 
gidity,the heavier the higher, but the flux of heavy nuclei 
is much less. Another models suggests that the cut-off ri- 
gidity of the main sources causes the knee. The present 
paper studies the spectrum and the anisotropy of cosmic 
rays in an isotropic diffuse model with explosive discrete 
sources in an infinite galaxyw 

2. The transport equation and its solution. 

In an isotropic diffusion model cosmic ray density, N(r, t, E ), 
obeys the following equation^ 

— — DV 2 N + —(NW) + BN *= Q. fi) 

dt dE 

where 0 is the diffusion coefficient, the third term is due 
to energy loss, ld=dE/dt, the fourth term is negitive source 
and Q is the source. Firstly, starts our argument from one 
discrete source and assume as follows: 

(l) Because the ages of the main point sources within lkpc 
(table. l) are about lO^y^j even for the highest energy 
proton its half lifetime( caused by Compton scattering, syn- 
chrotron radiation and bremstruhlung altogether) is several 
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decades to few hundred times long as its age. So we neglect 
the energy loss. 

(2) Assume the P-P interaction contributs to the negtive 
source term only* B^cn^ , where c is the speed of light* n 
the density of hydrogen* or^the inelastic cross-section of P- 
P interaction. From accelerator experiments^ the total cross 
-section of P-P is obtained among energy range 10-l*5xlCr GeV * 
and <^ e =0.8<^. Fitting these data* get the formula 

B(E)«3*10 -17 [35.3-2.2 ( lnE)+0.3< InE) 2 ] (1/s). (2) 
here n =1 has been taken. Table 1# List of SNR near the 
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(3) The charged cosmic 
rays propagate through 
resonant scattering with 
the turbulent bydromag- f8 <. 
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(3) 
is the 


NAME 

t(10 y) 

r(kpc) E (eV) 

CTB 72 

3.2 

0.7 

2x10 

Cyg. Loop 

3.5 

0.6 

10 

HB 21 

2.3 

0.8 

3x10 

CTB 1CT 

4.7 

0.9 

2x10 

CTB 13 

3.2 

0.6 

10 

HB g 

2.7 

0.8 

2x10 

S 149 

4.3 

0.7 

10 

Monoceros 

4.6 

0.6 

7x10 

Vela 

1.1 

0.4 

10 

Lupus 

3.8 

0.4 

4x10 

Loop 

3.0 

0.05 10 


12* e ABb 

where ¥ and 

mean intensity and energy 
density of the fundamental 
magnetic field* respectively. A is the fraction of turbulent 
in the mean field. Take Bs3^G, A s 0.01* ^=Q.3eV/cm 3 ,then for 
the proton of greater than lQGeV, 


D - 1.67 X 10‘ 5 E(cm 2 /s). 


I W 
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In general the coefficient has a form Do(E z * say E * but 
the recent data from HEAO-3 show that the index trends to 
rising and D°<E 0 ' 7 can account for the experimental data quite 
well. A compare shows that the diffusion coefficient adopted 
by our paper is close to that with the index 0.7 at high 
energies. 

(4) The i-th point source ejects particles transiently with 
a spectrum N 0 E* r at t=0. So that Q=N 0 E o(r-Y- ) J(t;) ^ 

Combining of all these and transforming N=J^e * we get 
solution of Eq.(l) _ r 


N i 


N o E 


r* 


(4irDt i ) 


13 exp[_ lB^~' B(E)t i 1 


(4) 


where t^ is the age of the point source* r^ the radial dis- 
tance with its origin at the source position. Mathematically, 
the formula( 4) makes sense for any value of r^ and t^, but 
considering the causality only such sources with (ri /t* )<c 
contribute to the observed flux. So we take account such 
sources only. 


3. The predict spectrum of proton. 

From Eq. ( 4 ) see that for very old/and distant pointsource 
its density N, has a very large Gaussian width, i.e.7N^=0, 
so we can treat them altogether as a background N^ . Therefor 
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we have 


N=N r + 

. & . v . 


supposing v is isotropy and neglecting higher than 
harmonic anisotropy of I, then N=4m/v, we have 

1=1 g+SL 47 A (6) 

Taking the parameters =1.5, N=4xl0 , Ig =2.5x10 m.s.sr .GeU ) 
and inserting these va- 
lues of r* ,t it of the ele— 10 <_ « 
ven supernovea in Eq.(6). vw° . o f 
we have the result(Fig.l), jp ° 3 

You can see that the pre- ^ 

diet results fits the « *. 

data of the proton quite % \ ^ 

well. In fact, because of « *' [b 

the cut-off energy of the 1= 10 
eleven supernovea, the 5 
value of Ig is almost uni- & 
quely deterraed by the 10 " 


second 


bw- 


data of the proton of 10 
-lO^eU. It is clear that 
the time scale of change 
of the proton intensity 
‘T>l/cn<£ e =lti 7 y,it agrees 
with long time data of 
cosmic rays. 

4. The anisotropy of 
cosmic ray. 

Under above mentioned 
case, the amplitude of 
first harmonic anisotropy 
has been obtained 

3r.N. • 

^ ( ^7 ) roi (7) 

where r t is the unit 
vector from the Earth to 
the i-th point source. 


there will be Compton— Ge- 
tting anisotropy, indepen- 
dent of energy but rele- 
vant to the spectrum 
index _ 

<£=(j+2)v/c (8) 

so the final resultant 
anisotropy, can be expre- 
ssed as 

(9) 
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Fig. 1 The differential spectrum of primary cosmic rays 
Curve 1 is the sum of intensities contributed by eleven superno- 
vae. The dashed line represents the background l Bt Curve 2 is 
the predicted differential spectrum of primary protons. Shaded 
area and sign Q are the data of protons, the rest are somewate 
contaminated by other nuclei. 
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Fig. 2. The expected amplitudes of the 
anisotropy caused by eleven SNR individ- 
ually. 
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The experimental data show that 
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N=47tl/c=1.2xlO" 3 E* ‘‘(1/m* ) for E$5xl0*GeV 
N=1 *1 E‘ 3 -‘(l/m*) for E>5xl0 GeU 

Inserting these values of r i ,t i , Nj,of the eleven supernovea 
individually and the value of N into Eq*(7), using the trial 
values N. o =1.32xl0 4 ' s , p=1.2(corresponding to the source emi- 
ssion 10 50 ergs in the particle form) , the anisotropies as a func- 
tion of energy were calculated. The results are given in Fig. 2. 

We find that each of these anisotropies is smaller than 
6x10 s below 10 ,4 eV, but varies as E 0-4 above 5x10' eV.So whatever 
the direction of each anisotropy is, the resultant anisotropy 
always has a E°' 4 variation. On hte other hand, the measured 
anisotropies are about 5.5x10*^ in the energy range 10"-10 l4 eV. 
Combination of these information together leads us to draw 
the conclusion that the Compton-Getting anisotropy is domi- 
nant, £=<£ =5.5x1 0**be low lCPeV.By substitutingtf=2.56 in Eq.(8), 
we get the velocity of 35km/s with respect to the cosmic ray 
background. 

5. Discussion and conclusion. 

The anisotropy mentioned above is in principle for proton. 
However, the anisotropy is dependent on the species of cosmic 
rays.Unf ortunately ,so far the identification of the species 
of cosmic ray in EAS is uncertain yet. Recently , a few of dis- 
cretej^-ray sources in the energy range 10 -10 ,6 e\/ has been 
detected.People have noted the possible effect of y”-rays on 
the anisotropy and the intensity of cosmic rays t9 i 

The conclusion to be drawn from this study may be summa- 
rized as follows: Taking some reasonable parameters, the model 
can account for the features of the proton spectrum and the 
approximate constancy of the cosmic ray intensity in a long 
period.lt also can inter prete the power law of the aniso- 
tropy above energy 5xlO ls eV ,and the Compton-Getting effect is 
responsible for the anisotropies in the energy range 10*'— IQ* 4 
eV . Furthermore, we got the streaming velocity of 35kra/s with 
respect to the cosmic ray background. 

Finally, I am grateful to Prof .A.W.Wolfendale for his helpful 
discussion and useful information. 
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